OBSERVATION OF THE CRAB NEBULA GAMMA-RAY EMISSION 
ABOVE 220 GeV BY THE CAT CHERENKOV IMAGING TELESCOPE 
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ABSTRACT 

The CAT imaging telescope, recently built on the site of the former solar plant "Themis" 
(French Pyrenees), observed 7-rays from the Crab nebula from October 1996 to March 1997. 
This steady source, often considered as the standard candle of very-high-energy 7-ray astron- 
omy, is used as a test-beam to probe the performances of the new telescope, particularly its 
energy threshold (220GeV at 20° zenith angle) and the stability of its response. Due to the fine- 
grain camera, an accurate analysis of the longitudinal profiles of shower images is performed, 
yielding the source position in two dimensions for each individual shower. 

INTRODUCTION 



With a reflector of only 17.7 m 2 , the CAT imaging telescope was designed to achieve a 
rather low threshold (~ 220 GeV) by taking full advantage of the rapidity of the Cherenkov 
signal with an almost isochronous mirror, fast phototubes and fast trigger and readout electron- 
ics. Moreover, its very-high-definition camera (546 pixels with an angular size of 0.12°) allows 
a powerful image analysis (Degrange and Le Bohec 1995) leading to a high rejection of hadron 
showers. A further advantage is that the source position on the sky can be reconstructed for 
each shower on the basis of the longitudinal distribution of light in its image, the accuracy 
being close to the pixel size. The position of a point source can then be determined to about 
1 arcmin with about 100 7-ray showers (Le Bohec 1996). 

The Crab nebula is the best-studied source of very-high-energy 7-rays (Lewis et al. 1993, 
Djannati-Ata'i 1995, Deckers et al. 1995) and its emission is found to be steady within the 



sensitivity of present instruments. It can thus be used as a test-beam to check the stability of 
the detector response and to characterize its energy threshold. Furthermore, it is important to 
validate the simulations of all instrumental effects (reflector optics, light catchers, phototubes 
and electronics) which are needed to determine the equivalent detection function of 

primary energy and zenith angle, as well as the accuracy in energy measurement, all quanti- 
ties necessary to measure the energy spectrum of the source. Such a validation is provided 
by comparing 7-ray event rates from the Crab nebula obtained experimentally with the cor- 
responding quantities predicted by simulations on the basis of the energy spectrum given by 
previous experiments, (Lewis et al. 1993, Djannati-Ata'i 1995). 

TRIGGER CONDITIONS 

The CAT triggering system uses a majority logic based on the 288 inner phototubes of 
the camera. The pulse heights of at least 4 phototubes were required to be greater than the 
average level of a 3-photoelectron signal, in at least one of 9 mutually overlapping sectors of the 
camera, each one corresponding to about 60°. With such conditions, random coincidences due 
to the sky noise are practically eliminated and the trigger rate, mainly due to hadronic showers 
and single muons, is ~ 15 Hz at moderate zenith angles (20° to 30°). These conditions do not 
correspond to the minimal threshold of the telescope, since the detector can be operated with 
trigger rates up to ~ 50 Hz. 

EVENT SELECTION AND BACKGROUND SUBTRACTION 

Data have been submitted to two different image analyses : 

Analysis based on the first and second moments of the light distribution 

The first method, based on the first and second moments of the light distribution in the 
image, has been successfully used by the Whipple group (Punch et al. 1991). Each shower 
image is characterized by its angular "length" L, "width" W, "distance" D (i.e. distance 
from the centroid of the image to the center of the field of view) and by its "size" S (i.e. the 
total integrated light content of the shower) (Fegan 1991). The cuts in these variables have 
been roughly adapted to the resolution offered by the CAT camera, but not optimized. The 
background due to hadronic showers as well as to single muons is severely reduced by the 
following requirements : 2 < L < 5 ; 0.7 < W < 1.5 ; D > 8.5 (in which angles are expressed 
in milliradians) and S > 30 photoelectrons. Then, the distribution of the orientation angle a 
(i.e. the angular deviation between the image principal direction and the line connecting the 
centroid of the image to the center of the field of view) shows a clear 7-ray signal at low values 
of a (a < 9°) above an almost flat background, which can be subtracted by using off-source 
data (figure la). 

Analysis based on a maximum-likelihood method 

The second analysis consists of a maximum-likelihood method (Degrange and Le Bohec 
1995) (Le Bohec 1996), based on an analytical model of 7-ray showers proposed by Hillas (Hillas 
1982). For a genuine 7-ray shower, this model provides the average distribution of Cherenkov 
light in the focal plane as a function of several parameters : the azimuth angle </> of the image 
axis, the 7-ray energy £" 7 , the distance R between the shower axis and the telescope (or impact 
parameter), and finally the two angular coordinates of the source. In this method, the source 
position is thus left as a free parameter which can be reconstructed on a shower per shower 
basis, due to its influence on the longitudinal light distribution. The analytical model is used 
to define a x 2 -like function of <fi, E 7 , R and source position. Fluctuations in the pixel signals 



are estimated from Monte-Carlo studies of 7-ray-induced showers by requiring that the x 2 
probability distribution be approximately uniform. For each shower observed, this function is 
minimized with respect to the parameters: the fitted % 2 provides gamma-hadron discrimination; 
the 7-ray energy determination automatically takes account of the shower position with respect 
to the telescope which is obtained from the same fit, which further yields the source position 
in 2 dimensions. For those events with a \ 2 probability greater than 0.35, the distribution 
of the orientation angle a is shown in figure lb. Since the method discriminates between the 
top and the bottom of showers, this angle is here plotted from 0° to 180°. The 7-ray signal 
clearly shows up at low values of a (a < 9°) whereas the almost flat hadronic or muonic 
background can easily be subtracted by using off-source data. The 2-dimensional distribution 
of the reconstructed source position is shown in figure 2 in celestial coordinates centered on 
the Crab nebula, after background subtraction, showing that a single imaging telescope with a 
very-high-definition camera has a real source localization capability. 

GAMMA-RAY EVENT RATES 

The simulation is based on a two-step procedure, one concerning the development of 
electromagnetic air showers (Kertzman and Sembroski 1994) and the other the detection of 
Cherenkov photons and subsequent instrumental effects (Le Bohec 1996). In the present study 
whose goal is to validate our knowledge of the detector response, the Crab nebula spectrum 
quoted in Lewis et al. 1993 was used as input to simulate 7-ray-induced showers above 50 GeV 
with an integral spectral index of 1.7. Losses in 7-ray events due to the cuts described above 
can easily be modeled in the moment-based analysis, whereas the second method requires a 
more accurate treatment of the night-sky background; we thus restrict ourselves here to the 
first method. For those observations with a zenith angle smaller than 30°, the raw 7-ray rate 
(i.e. after cuts) from the Crab Nebula is found to be 2.0 ±0.3 7-ray per minute. Correcting for 
losses due to cuts, the total rate is found to be 5.1 ± 0.7 7-ray per minute. This corrected rate 
characterizes the telescope response and particularly its energy threshold. Since the definition of 
the threshold of an Atmospheric Cherenkov Telescope varies among authors, the corrected 7-ray 
event rates from a steady source provide an objective way of comparing different instruments. 
The Monte-Carlo expectations for the corrected event rates (6.4 7-ray per minute at 20° zenith 
angle and 4.77-ray per minute at 30°) are found in good agreement with the experimental value, 
showing that the instrumental simulation provides a reasonable description of the apparatus. 
If the threshold is defined as the energy of maximal 7-ray-event rate from the Crab nebula, 
the corresponding value for CAT with the triggering conditions quoted above is 220 GeV at 
20° zenith angle. Finally, figure 3 shows that 7-ray event rates obtained in similar conditions 
of zenith angle but at different times are compatible with each other within statistical errors; 
the detector response is thus found to be stable, a necessary condition for the study of very 
variable sources such as the Active Galactic Nuclei Markarian 421 and Markarian 501. 

CONCLUSION 

From the preceding study, we conclude that the CAT telescope, which started operation in 
October 1996, achieves a low energy threshold despite its moderate size, as shown by the 7-ray 
event rates from the Crab nebula. Its response is found to be stable in time. Furthermore, its 
very-high-definition camera provides a real source localization capability, which will be essential 
in the study of poorly localized sources. 
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Fig. 1: Orientation angle a for 8 hours of Crab observation : (a) moment based analysis and 
(b) maximum-likelihood method. 
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Fig. 2: Reconstructed source position after background subtraction in celestial coordinates for 
the Crab nebula. 
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Fig. 3: 7-ray event rates for 90 minute periods in similar conditions of zenith angle (20° - 30° , 
December 1996 to Febuary 1997). Horizontal scale is arbitrary. 



